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1. Introduction 

Exact solutions of the Einstein equations play a very important role in cosmology, because 
they permit to study in a convenient way the qualitative and quantitative features of the 
behavior of the universe as a whole. During the last decade the cosmological models with 
scalar fields have acquired a great popularity. It is worth mentioning various scenarios of 
the inflationary expansion of the early universe Q or the the quintessence models of the 
dark energy Q responsible for the phenomenon of cosmic acceleration Q . Nevertheless the 
number of known exact solutions for cosmological models based on scalar fields is rather 
limited. One of such models is the flat Friedmann universe filled with a minimally coupled 
scalar field with exponential potential. A particular solution for this model was known 
since the eighties and was studied in detail 0, |5|, ^, ^ |9|. This solution describes a 
power-law expansion of the universe. More recently, the general solution of the Einstein 



equations for this model was constructed [^-|18]. This general solution was used for the 
description of such effects as transient acceleration, and for the analysis of some models 
related to strings and branes. 

Notwithstanding this activity, a good description of the general solution, emphasizing 
its difference from the "old" particular solutions |Q, ^, 0, ^, ^, is not available to our 
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knowledge. In addition, all the considerations of the general solution for the scalar field 
with an exponential potential involved a real scalar fields with real potentials. However, 
there are attempts to consider cosmological models with complex scalar fields and complex 
potentials. Such models do not contradict to common sense, if the observable (first of 
all, geometrical) characteristics are real. In our preceding papers inspired by the 
development of the PT-symmetric quantum theory |2l[], we have developed models with 
complex potentials and have shown that they are rather convenient for the description of the 



so called phantom cosmology [23| , including such an enigmatic phenomenon as phantom 



divide line crossing [l24 |. 



In the present paper we provide a comprehensible description of the general cosmo- 
logical solution with an exponential potential; then we explore the general solution for the 
phantom field. Finally, we describe the general exact cosmological solution for the case of 
a piecewise exponential potential. 

The paper is organized as follows: in Sec. 2 we describe the general cosmological 
solution with an exponential potential in contrast to the "old" particular solution; in Sec. 3 
we present the phantom version of this solution; in Sec. 4 we construct the general solution 
for a piecewise exponential potential with cusps; the last section is devoted to Conclusions. 

2. The general solution with an exponential potential versus the "old" 
particular solution 

We shall study the flat Friedmann cosmological model described by the metric 

ds"^ = dt^ - a^{t)df, 

where a{t) is the cosmological radius of the universe. The dynamics of the cosmological 
evolution is characterized by the Hubble variable 

h= -, 
a 



which satisfies the Friedmann equation 



h^ = e, (2.1) 



where £ is the energy density of the matter populating the universe and "dot" means the 
derivative with respect to the cosmic time parameter t. 

We now basically follow the approach presented in Ref. ||l6[| , adapting it to our pur- 
poses. We shall consider the flat Friedmann universe (^) fllled with the minimally coupled 
scalar fleld with the potential 

V{(t>) = Voe^t (2.2) 
Now the Friedmann equation ( |2.1| ) has the form 

-,=^ + Voe't (2.3) 

while the Klein-Gordon equation is 

(?i + 3-0 + AFoe^* = 0. (2.4) 
a 
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2.1 The particular exact solution 

A particular exact solution |5|, ^, ^, |8|, describing the power-law expansion of the 
universe can be found as follows: Suppose that the scalar field has a time dependence 

(t>{t) = ^olnt + ^i, (2.5) 

while the cosmological radius behaves as 

a = aot''. (2.6) 
In this case the Klein-Gordon ( |2.4| ) equation acquires the form 

_^+W» + ^=„. (2.7) 



It follows immediately that 



Then 



2 

A' 



(2.8) 



2(3A: - 1) = AVoe^-^i. (2.9) 
The Priedmann equation ( |2.3| ) gives now 

k"^ = ^ + Voe^'f". (2.10) 

Combining Eqs. (U) and ( ^lo|) we obtain 

k = -^. (2.11) 

Substituting ( p. 11 ) into the relation ( 2.10| ) and requiring positivity of the coefficient Vq, we 
obtain the following restriction on A: 

|A| < 3\/2. (2.12) 

(Note that the limiting case |A| = 3V2, would imply k = 1/3, i.e. the law of expansion 
of the universe filled with stiff matter or massless scalar field, which in turn, means that 
Vo = 0). We can find also the value of the constant 0i: 

Here, let us note that for the case of an imaginary A the particular solution still exists 
and is purely imaginary [19|. Here, both the potential and kinetic energy are real, but the 
latter is negative and, hence, we encounter the phantom kind of matter. 

It is convenient to write down again the explicit form of the exact particular solution: 
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h{t) = ^. (2.15) 

The important feature of this solution is the rigid relation between the (/>(t) and its time 
derivative (j){t): 

2 / 2 \ 1 /2(18-A2)\ , , 

Thus, the particular exact solution corresponds to a particular choice of the initial condi- 
tions: at any moment of time fixing the value of (p, we automatically fix the value of cj). 
It means that on the phase space of the problem under consideration the particular exact 
solution is described by a unique trajectory. It is convenient to introduce a new variable: 



^> = VVoe^'t'. (2.17) 
Choosing the phase space variables as <I> and (p we can calculate their relation: 

6 2A 



^ \/36 - 2A2 



const. (2.18) 



Thus, the trajectory on our phase space is a ray, beginning at the infinity and ending at 
the coordinate origin. Correspondingly, the universe begins its evolution from a Big Bang 
singularity and then undergoes an infinite expansion (seeEq. ( ^l5|) . 

2.2 The general solution 

Now, let us turn to the construction of the general exact solution. It is convenient to 
introduce new variables, u and v such that 

= 6^^+", (2.19) 

(l) = A{v-u), (2.20) 
where A is a coefficient to be defined. Now the Friedmann equation ( |2.3D is 

-(«2 + {,2 + 2uv) = —{u^ + {;2 - 2uv) + Foe^^^^""). (2.21) 
9 2 



The Klein-Gordon equation ( |2.4D has now the form 

A{v -u) + A{v^ - ii^) + AFoe^^^^-''^ = 0. (2.22) 
Choosing the coefficient A as 

A=^, (2.23) 
we give a "light-cone form" to Eq. ( 2.21| ): 

vu=^Voe'T^^''-''). (2.24) 

Now we want to simplify further this equation, choosing a new time parameter r. Rewriting 
Eq. (1^ as 

v'u'f^ = ^^oe'?^^^-") (2.25) 
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where prime denotes the derivative with respect to r. It is convenient to choose this new 
time parameter such that 



we come to 



3 ^ 3 / Hv-n) 



V u = 1. 



Correspondingly, the Klein-Gordon equation (2.22) is now 

// // I V2A , i\2 , / /2 i2\ , 2-v/2 . 

V -u H —{v - u) +{v -u )H — A = 0. 

fa 6 

Substituting u' = 1/v' from Eq. ( |2.27D into Eq. (|2.2g| ) one arrives to 



v"+\l + 



V2X 



Introducing a new variable 



/2 I 

V + 



X = V 



V2X 



1=0. 



we rewrite the Klein-Gordon equation ( p. 29 ) in the Riccati form: 



x' + 1 + 



V2X 



V2X 



1=0. 



(2.26) 
(2.27) 

(2.28) 

(2.29) 
(2.30) 

(2.31) 



Making the substitution 



1 



y 



1 + 4^) y' 



we obtain the second-order linear differential equation 



(2.32) 



(2.33) 



We should consider separately two cases : the "hyperbolic" one, when the constant A 
satisfies the condition ( 2.12 ) and "trigonometric" one, when the condition ( p.l2| ) is not 
satisfied. 

2.2.1 The hyperbolic case 



The solution of Eq. ( p.33|) in this case is 

y(r) = fie"" + Ce"''", 

where 

Then, substituting ( 2.34| ) into Eqs. ( |2.30| ) and ( |2.32| ) after an elementary integration we 
find 

■ln{Be^^ + Ce-^^) + vo. (2.36) 



(2.34) 



(2.35) 



1 



1 + 



6 ) 
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Using the relation ( 2.27] ) we can analogously find 

1 



u 



1-^ 

Now, remembering the formulae ( |2.2C1| ) and ( p. 23 ), we can write 

1 



(2.37) 



V2 

^ I (^^0 - uo) + 



1 



V2\ \ 

6 



1 + 



Ce- 



(2.38) 



We have to consider separately two cases: C = and C ^ 0. If C = the solution 



( 2.38 ) can be rewritten as 



9av 



where (p2 is an arbitrary constant. Substituting Eq.( p.39| ) into Eq. (2.26), connecting the 
time parameter r with the cosmic time t we obtain the following equation 



dT 3 /- 



A(/>2 



9k 



(2.40) 



Integrating the equation ( 2.4C| ) we find the parameter r as a function of t and substituting 
it into Eq. ( p9D we come to the particular exact solution ( |2.14| ). The dependence on the 
constant (/)2 disappears. Thus, the case C = coincides with the particular exact solution 
described in the subsection 2.1. One can easily show that the case 5 = corresponds to 
the particular solution, described in the preceding subsection, but it refers to a universe 
contracting from a state with an infinite radius to the Big Crunch singularity. 
In the case when C / 0, B / the solution (|]3|) can be rewritten as 



2At 
"9k 



+ 



V2 



3(1 + 4^) 



ln(l + ^e-2- 



3 1 



^mfi-^e-- 



(2.41) 



where (p^ is a constant. As a matter of fact, here there are two families of solutions. We 
can rescale the factor ^, that is equivalent to a shift of the variable r. Let us consider first 
the case C / B = 1. Then we come to a more simple expression: 



(t> = (t>3 



2Xt 



+ 



V2 



3(1 + ^) 



ln(l + e 



-2kt\ 



V2 



ln(l - e'^""^). (2.42) 



3 1 



V2X\ 
6 ) 



While the particular solution ( 2.14| ) describes the unique trajectory in the two-dimensional 
phase space of variables (<^, </>) which can be represented by Eq. ( p.l6| ) or by Eq. ( 2.18| ) the 
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general solution ( 2.42| ) describes a family of solutions, parameterized by the constant 4)3. 
Indeed, in this case a rigid relation between (f> and (p is absent. If one fixes the value of <j) 
then the corresponding value of (f> is given hy (j) = 4>' -f. As it follows from Eq. ( p. 26 ) fixing 
the value of (j) determines the value of f. However, changing the value of the constant ^3 
while preserving the value of (j) is equivalent to shifting the value of r, and, hence, changing 
the value of (j)' and (p. One can represent it in a different manner. Let us note that the 
value of r in Eq. ( p. 42 ) is changed in the interval < r < 00. We can now calculate the 
relation between (j) and the phase space variable $ analogous to ( p. 18 ). One obtains 



A 
3k 



\/2e 



-2ftr 



k{1 



2 + 



\/2A _ 



2kt 



(2.43) 



Note, that the value of both (j)' and (j) is negative during the cosmological evolution. 

In contrast to the formula ( p. 18 ) the formula ( |2.43 ) does not describe a ray. The 
corresponding curve begins at the point <I>(t = 0) = +oo,0(r = 0) = —00 and ends at 
r = +00 at the origin of the coordinates. It is easy to see from Eq. (gl|) that the 
value of ^ at the moment r = tends to — c«. That means that at the beginning of the 
cosmological evolution the kinetic term dominates the potential term. It is well-known 
that in this case (the case of the massless scalar field or, equivalently, the stiff matter) the 
Hubble parameter behaves as 

/i(t)^l t^O. (2.44) 



The direct evaluation obtained by the substitution of the exact solution ( 2.42 ) into the 
Friedmann equation (|2.3| ) and using the relation ( 2.26| ) in the limiting case t — )■ confirms 
the behaviour ( p. 44 ). At the end of the evolution r = +00, the relation ( p. 431) coincide 
with that of the ray ( 2.18| ). Thus, the curves, corresponding to the general solution ( 2.42| ) 
(excluding the particular solution ( |2.14| )) begins at the Big Bang singularity, where the 
kinetic term dominates the potential one, which in the phase plane ($, cf) corresponds to 
the asymptote coinciding with the negative semi-axis $ = 0. Then, the curves ( 2.431 ) are 
located under the ray ( 2.18| ) and conclude their evolution at the coordinate origin, which 
corresponds to an eternal expansion. 

Now we can consider the second family of solutions, which can be obtained by putting 
C/B = —1. The time parameter again runs between and +00. However, the expression 
for the scalar field is now 



2At 



+ 



V2 



ln(l 



-2kt 



^/2 



3 1 



V2\ \ 
6 ) 



ln(l + e 



-2kt\ 



(2.45) 



This solution also describes an evolution from the Big Bang to an infinite expansion, 
however, now at the beginning of the evolution the scalar field (p — )■ —00 and the potential 
energy is equal to zero. During the evolution, the scalar field is growing, then at some 
moment its time derivative changes the sign and it becomes decreasing. At the end of the 
evolution (an infinite expansion) when r — )• cxd the scalar field again tends to —00 and both 
the kinetic and potential energy vanish. 
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The particular exact solution plays a role of the separatrix between these two families 
of solutions. One can ask how to understand to which exactly family belong a certain 
trajectory being given the values of the scalar field (j) and its time derivative (with respect 
to the cosmic time t) at some moment of time. The answer is simple. It is enough to 
calculate the derivative of the scalar field with respect to the parametric time r, which is 
(j)' = ^— If its value is less than — then the trajectory belongs to the first family. 

In the opposite case it belongs to the second one. In Figure 1, we sketch the trajectory 
describing the particular exact solution and two trajectories, representing two families of 
the general solution for the hyperbolic case (for definiteness, we have chosen for all figures 
the subcase A > 0). 




Figure 1: The phase space diagram $,7r, where tt stays for the velocity tt = </>. The straight line 
represents the exact particular solution for the hyperbolic case, while two curved lines represents 
solutions, belonging to two families of general hyperbolic solutions described in this subsection. 



2.2.2 The trigonometric case 

Now we consider the case when > 18. In this case the solution of Eq. ( 2.33| ) can be 
chosen as sinwr or cosa;r. However, in contrast with the hyperbolic case, these two choices 
do not imply the real difference between cosmological evolutions, because the transition 
from one to another corresponds to the shift of the parametric time interval. Thus, we 
shall choose 

y(r) = DsincjT, (2.46) 

where 

- = /f^' (2-47) 
and D and ao are real constants. Then using the definitions of ( |2^ and ( ^^301 ) we get 

v = ^ In sin wr + , (2.48) 

1 _i_ v2A 
^ 6 



and similarly 



—= Incoswr + UQ. (2.49) 

v2A _ 1 
6 ^ 
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Finally, the general solution for this trigonometric case looks like 




InsmwrH — -= Incoswr , (2.50) 

V 2A _ 1 / 
6 / 

where (\)^^ is a constant. We note that the expression describes an expanding universe if the 
parametric time r is changing inside the interval 

0<a;r<|. (2.51) 

The relation between the parametric time r and the cosmic time t is 



2 

We shall present also 



3 /^gA04/2 (g-j^^^) 1+4^ (pQg^^) 4^-1^ (2.52) 



2x/2(4A 



cos 2CJT — 1 



(2.53) 
3a; sm 2ujt 

Using the formulae (IZSOl) , and (|I5^ ) one can find that when r runs from to 7r/2a;, 



the cosmic time runs from to oo. At the beginning of the evolution the velocity ^ is 
infinite, while the potential energy is zero. The direct calculation shows that the Hubble 
parameter behaves at the beginning of the evolution as = ^ as it should be for universe, 
filled with a massless scalar field arising from the Big Bang singularity. At the end of 
the evolution the potential energy again tends to zero, while the velocity (\) being negative 
also tends to zero. (Note that the time derivative with respect to the parametric time r 
diverges when lot tt/I and one can think that here we again encounter the singularity as 



was mentioned in |16|. However, it is not the case, because the expression f tends to zero 
more rapidly and thus, the kinetic energy of the scalar field tends to zero). This situaiton 
corresponds to an eternal expansion. From Eq. ( |2.53| ) it follows that the velocity changes 
sign at the moment 

1 3^2 

Ti = — arccos — - — . (2-54) 
luj A 

In Figure 2 we represent a typical trajectory for the trigonometric case. 
2.2.3 The case A = \/l8 

In this case the equation for the variable v is the following : 

v" + 2v''^ = 0, (2.55) 

whose solution is 

v = ^\iiT + vo. (2.56) 

Correspondingly 

u = t'^ + U0 (2.57) 
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Figure 2: The phase space diagram $, tt, where tt stays for the velocity tt = 0. The curved 
line represents one of trajectories, belonging to the family, given by the general solution for the 
trigonometric case. 



and 



<P5 + ^iy^li^r-T']. (2.58) 



The relation between r and t is now 



f = ^y%e^^^^/^V^e-^\ (2.59) 

It is easy to see that r runs from to oo and the cosmic time t has the same range. Then 

V5 _ 2V2r 

^ 6r 3 ^ ^ 

It is easy to see that at the beginning and at the end of the evolution the field (p tends to 
— oo and hence the potential energy vanishes. Then, at the beginning of the evolution the 
(f) is infinite and that means that the universe is born from the Big Bang singularity, while 
at the end of the evolution the negative time derivative of the scalar field (p vanishes and 
the universe is expanding infinitely. 

It is curious to note that the particular solution described in the subsection 2.1 exists 
in the case = 18. In this case, as it follows from Eqs. (PD-(^1ID, 



Vo = (2.61) 

and this is simply the case of the massless scalar field. This particular case cannot be 
extracted from the general solution ( 2.5^ ) because the very existence of this general solution 



is based on the introduction of the parametric time r which is defined by means of Eq. 



( 2.26 ) which is senseless when Vq = 0. Thus, the particular solution for the hyperbolic 
case is contained in the general formula ( |2.38| ). In the case = 18 the particular solution 
exists, but it should be treated separately, while the trigonometric case does not have a 
particular solution. 
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3. The exact solution for a phantom scalar field with an exponential po- 
tential 



Considering the phantom scalar field with a negative kinetic term, we shall have the fol- 
lowing Friedmann 

^^--f + l'oe^^ (3.1) 

and Klein- Gordon 

4> + 3-<P- XVoe^'f' = (3.2) 
a 

equations. We shall first present the particular solution for this case. 
3.1 The particular exact solution for the phantom case 

As in the subsection 2.1, we shall look for the solution for the phantom scalar field, which 
depends logarithmically on the cosmic time t, while the scale factor undergoes a power-law 
expansion (or contraction). If we would like to consider an expanding universe, then the 
solution will be 

„<)._^M-.) + iln(^). ,3.3) 

where t is running from — oo to 0. The Hubble parameter is now 

hit) = -^^. (3.4) 



Thus, the formulae (3.3) and (3.4) describe a cosmological evolution which begins at t = 
— oo with infinitely small cosmological radius and ends at t = 0, encountering a Big Rip 
singularity. However, another particular solution for the phantom case does exist. Here 
the formula for the Hubble factor is the same (p.4|), but in the formula for the scalar field 
( |3.3| ) Int instead of ln(— t) enters. This solution describes the contraction of the universe, 
which begins at the moment t = in the "anti-Big Rip" singularity characterized by an 
infinite radius and infinite negative Hubble parameter, and ends a.tt = +oo with an endless 
contraction. 

3.2 The general exact solution for the phantom case 

As in the subsection 2.2 we shall introduce the variables v and u (see Eqs. (2.19),( 2.2[1| ), 



( 2.23D ). However, because of the negative sign of the kinetic term in the right-hand side of 



Eq. {^) we shall obtain instead of ( |2.24| ) 



y^ + u'^ = ^V^e^'l>. (3.5) 



It is convenient now to introduce a complex variable 



Now, Eq. (|3^) looks like 



z = —;=(v + iu). (3.6) 
V2 



z'z = Ivoe^''', (3.7) 
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where "bar" stands for the complex conjugation. Introducing now the time parameter r 
hke in Eq. (2.26) we obtain 

(3.8) 



z'z' = 1. 



Rewriting the Klein-Gordon equation (3.2) and taking into account the relation ( |3.8| ) we 
come to 



1 + 



V2Xi 
6 



V2\ 
~6~ 



0. 



If 



Introducing the function / such that 

^/ 

a f ' 

V2{l-i) ( y/2\i 



where 



a 



1 + 



6 



The auxiliary function / satisfies the following equation 



f" - f^'f = 0, 



where 



A2 

1 + 



The general solution of Eq. ( 3. 121 ) is 

/(r) = Fe~^^ + Ge- 

Hence, 



z = - Inf. 

a 



Now, we can find the expression for the scalar field 



, , 2A . . 2\/2 
v) = TT^^^I/I ~ ^T^arg/ + const. 



(3.9) 

(3.10) 
(3.11) 

(3.12) 
(3.13) 

(3.14) 
(3.15) 

(3.16) 



Now, substituting the expression for / ( 3.14| ) into the expression for z ( p.l5| ) and substi- 
tuting the latter into the condition (^]^) we come to the consistency equation 



FG + GF = 0. 



(3.17) 



These equation can be satisfied if one of the coefficients is equal to zero or if the difference 
of their phases is equal to 7r/2. 

First, we consider the case G = 0. Then 



(3.18) 



Substituting this expression into Eq. (|2^ ) we can find r as a function of the cosmic time 



A2 



ln(-t) + ^ + ln^ ,~ 

2 \ OK 
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Substitung this expression into the solution ( ^.18 ) we come to the particular exact solution 



( |3.3| ) as it should be and in complete analogy with the hyperbolic case. It is easy to see 
from Eq. ( p. 19 ) that when the cosmic time runs from — oo to 0, the parameter r runs 



from — oo to +oo, while the cosmic time t runs from t = — oo to t = 0. In contrast to the 



hyperbolic case considered in subsubsection 2.2.1 (cf. Eq. ( 2.38| )) we can here put also 
F = 0, while G ^ 0. In this case we reproduce the second particular solution, describing 
an infinite contraction of the universe, which begins at the anti-Big Rip singularity. 

Now, we consider the case when both the constants F and G are different from zero. 
We can choose one of these constants, say F equal to 1 while G = i. 



Ho + 7^1ncosh2Kr - ^j^arctan(e~^^^). (3.20) 



2A , , 2^2 
3^ 

Here, the time parameter r runs from — oo to +oo. When r — >• — oo the field (j) — t- +oo, 
while when r — )• +oo the field (j) ~^ +c« again. Thus, the beginning and the end of the 
cosmological evolution are characterized by a positive infinite value of the scalar field, and, 
hence by the positive infinite value of the potential V. It is useful to write down also the 
explicit expression for the Hubble parameter: 

Kr) = -~ ^ . 3.21 

y cosh2Kr j 

It is easy to see that at r — )• — oo the Hubble parameter h — t- — oo at r — )■ +oo the Hubble 
parameter h — )• +oo and at r = ^^arcsinh-^^ the Hubble parameter changes the sign 
and, hence, the universe passes through the point of minimal contraction. The complete 
cosmological evolution involves a finite period of the cosmic time t. One can say that these 
finite time evolutions represent the phantom counterpart of the well-known cosmological 
evolutions, which begin in the Big Bang singularity, reach the point of maximal expansion 
and then have a stage of contraction culminating in the Big Crunch singularity. These 
evolutions can be observed not only in the closed Friedmann models, but also in the fiat 
Friedmann models with the standard scalar field with a negative potential. (As far as we 
could understand this solution beginning from "anti-Big-Rip" singularity and ending in the 
Big Rip singularity, passing through the point of minimal contraction was not considered 
in paper [17| devoted to the study of phantom solutions with exponential potentials). 



In Figure 3 we represent two particular exact solution and a typical example of general 
solution for the phantom case. 

3.3 The complexification of the scalar field and of its potential 

In principle, the phantom solutions can be obtained from the hyperbolic non-phantom so- 
lutions by some kind of analytic continuation. Let us consider the case when the parameter 
A = iA (where A is real) is imaginary. One can see that this case is a ramification of the 
hyperbolic case, because the parameter k = VlS -|- A^ ( ^.35 ) is well defined. However, now 



also the solution of Klein-Gordon equation should be purely imaginary. Thus, taking Eq. 



( 2.33 ) with K, = \/l8 + A^ we obtain the general equation of the hyperbolic type, which can 



be written as was don in Eq. (2.38) with arbitrary coefficients B and C. Now, choosing 
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Figure 3: The phase space diagram $, tt, where tt stays for the velocity tt = (j>. Two straight 
Hne trajectories describe two particular solutions, corresponding to the expanding and contracting 
universes. The curved line represents a trajectory, belonging to the family of those evolving from the 
anti-Big Rip singularity to the Big Rip singularity, passing through the point of minimal contraction 
of the universe. 

these coefficients as i? = 1 and C = i, we come to purely imaginary solution, which behaves 
just like our phantom solution ( [j.20D . multiplied by i. Substituting this solution into the 
Friedmann equation we obtain an equation of the evolution of the universe, filled with a 
phantom scalar field. Naturally, all the observables appear to be real, and when necessary 
also positive (cf. []T9|| ). 

4. General exact solutions for piecewise exponential potentials 

In this section we analyze the generalization of the previous results for the case of the 
potentials, represented by piecewise functions, where all the pieces are exponential. We 
limit ourselves by continuous (but not smooth) potentials. Such potentials can be written 
down as 

+9{^-(l)i)ei(j)2-<P)Vie^'''' + --- 

+e{cl) - <PN)VNe^^t (4.1) 
where 9 is the Heaviside function and where 

-oo < (1)1 < 4>2 <■■■< <f>N < oo, (4.2) 

and 

y^gAfe0fc+i ^ Vfe+ie-^'-'+i'^'^+S < A; < iV - 1, (4.3) 

The general exact solutions exist also in this case and can be realized by means of the 
matching across the values of scalar field, where the exponent A changes its value. This 
matching involves not only the value of the scalar field, but also the continuity of (p and, 
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hence, (p' . In the opposite case, one would encounter a jump in the value of the Hubble 
parameter, following from the Priedmann equation, which seems unphysical. The simplest 
potential with cusp of this type is V{(j)) = Vqc'*''''^!, which can be obtained from our general 
potential by choosing = 1, 0i = 0, Vi = Vq and Aq = — Ai = —A. 

To get an idea of new possibilities let us consider, for example, the matching across the 
value (f) = (j)i and focus on an expanding universe. If Aq satisfies the relation Ag > 18 (the 
trigonometric case) two situations are possible: the derivative (j)' can be positive (an initial 
part of the evolution) and it can be negative (the final part of the evolution) as follows 
from Eq. ( 2.531) . If the derivative is negative, that means that we cross the value (j) = 
from the right ( from the side of larger values of (p) and then the value of is decreasing 
indefinitely. 

If this derivative is positive the field crosses the point (f) = (f)i and enters the range of 
values between (pi and (p2, where the exponent A is equal to Ai. If Ai is also trigonometric, 
then the universe continues its expanding evolution, with the scalar field given by Eq. 
(2.50) with A = Ai and the corresponding initial conditions. Then, again two situations 
are possible - the field continues growing arriving to the value (j) = 4'2i where another 
change of regime occurs, or it can arrive at some moment to the value (p = (p.^ < (p2 while 
(p' = 0. After that the field begins decreasing arriving at the value (p = (pi with a negative 
derivative (p' and enters into the range of values with (p < (pi, decreasing indefinitely until 
— oo, which corresponds to an eternal expansion with the vanishing value of the Hubble 
parameter. 

If the value of A^ < 18 (hyperbolic case) and if the value of (f>' is less than — |^ then 
the solution in the range 0i < < 02 is the hyperbolic solution of the first kind (described 
in the subsubsection 2.2.1) and the field enters into the region (p < (pi from the right, 
continuing to decrease indefinitely. If the value of p' at the point, where (p = pi is negative 
but larger than — then the solution in the region </)i < </> < </>2 is the hyperbolic solution 
of the second kind and again the scalar field after the crossing the value (p = (pi decreases 
indefinitely. If instead the value of (p' at the point of transition is positive the fields enters 
into the region 01 < </> < 02 increasing and following the hyperbolic solution of the second 
kind. Then, again, two situations are possible: the field p growing can achieve the value (p2 
with the subsequent change of the regime or it can begin decreasing after achieving some 
maximal value (pmax < 02 , entering back into the region < 0i. 

It is important to pay a special attention to the situation when the field arrives at 
the value = 0i with 0' = 0. As follows from Eq. (2.53) that means that the field has 
reached its maximal value and begin decreasing. Thus, in this case the change of regime 
does not occur. In our preceding papers |25, ^ (see also ||2^), in was shown that the 
phantom divide line crossing effect, provoked by self-conversion of a non-phantom scalar 
field into a phantom scalar field (or vice versa) can take place in the models with a unique 
minimally coupled scalar field. To realize this effect it is necessary to consider the cusped 
potentials and some particular initial conditions. In other words, it is necessary to approach 
the cusp with a vanishing time derivative of the scalar field 26 1. We see that in the 
case of piecewise exponential potentials such an effect is impossible, in spite of the presence 
of cusps. As a matter of fact the cusps, permitting such an exotic phenomenon should be 
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non-analytic |25, while the cusps, which we treat now are still too smooth. 

Concluding this section, we remark that choosing some simultaneous values of (p and 
(f), or, equivalently (j) and (j)' and having the piecewise potential (^]^) one is able to recon- 
struct the whole past and future evolution of the universe, using the general exact solutions, 
described in section 2. One can describe also the evolution of the phantom field in the piece- 
wise potential using the results of section 3. However, the effects of (de)-phantomization 
do not occur. As was mentioned above, the effect of (de)-phantomization requires the 
presence of strong non-analyticity of potential. Besides, the field should approach the cusp 
with zero velocity and cross it. In the case considered in this section instead of crossing we 
have a reverse motion. 



5. Conclusions 

As far as we know the only general exact cosmological solution in the presence of a scalar 
field, constructed explicitly, is the solution for a flat Friedmann universe filled with a 
minimally coupled complex field with exponential potential. While a particular exact 
solution for this case was known from 1985 0, the general solution was constructed only 



in 1998 [10|. In the present paper, we have studied accurately this solution, correcting some 
imprecisions, encountered in the literature. In addition, we have generalized this solution 
for a piecewise exponential potential, which is continuous, but not smooth (with cusps). 
In spite of presence of cusps, in such a model the effect of (de)-phantomization does not 
exist p^, 26 1 which requires non-analyticity of the potential. 
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